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Abstract 

We show that the short-time critical exponent 6 related to the critical initial slip in a stochastic 
model can be determined by the time correlation of the order parameter. In our procedure it 
suffices to start with an uncorrelated state with zero order parameter instead of departing, as 
usually done, from an initial state with a nonzero order parameter. The proof uses the group of 
symmetry operations related to the Markovian dynamics. Our scheme is extended to cover models 
with absorbing states. 

PACS numbers: 05.70.Ln, 05.50.+q, 64.60.Ht 
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I. INTRODUCTION 



The universal behavior occurring in the first steps of a Monte Carlo simulation, the short- 
time dynamics, has been amply investigated in the last years. According 
to renormalization group arguments Q the early time behavior of the order parameter (the 
magnetization, for example, in the case of the Ising model) follows a power law with a 
critical universal exponent 9. The numerical calculation of the exponent 9 is performed by 
placing the system at the critical point and departing from a configuration where the order 
parameter mo is very small. 

The purpose of this article is to show that it is possible to determine the exponent 
9 by starting with a configuration in which the order parameter is identically zero. In 
our approach we do not measure the order parameter itself, which is zero, but its time 
coorelation. If M(t) denotes de instantaneous order parameter we show that the quantity 

Q(t) = (M(t)M(O)) (1) 

follows a power law 

Q(t) ~ t e (2) 
where the initial configuration is uncorrelated and such that 

(M(0)> = (3) 

This result is general and can be applied to any lattice system described by a Markovian 
process and such that the transition probability is invariant under a given group of symme- 
try operations. In this procedure, which was previously deduced for models with up-down 
symmetry ^]. Here we generalize the scheme to indued models with other types of sym- 
metries. In particular we apply the method to systems with antiferromagnetic ordering 
and to the Baxter- Wu model in a triangular lattice [9|]. Employing this procedure it was 



recently calculated by numerical simulations lOj], with an excellent precision, the exponent 



9 associated to the Baxter- Wu model with triplet interactions in a triangular lattice. 

We consider also the short-time behavior of intrinsic irreversible models with absorbing 
states, such as the contact process. In this case, we demonstrate here that it suffices to depart 
from a configuration with just one occupied site instead of starting from a configuration 
where the density of particles is finite and small. Therefore, the short-time behavior of the 



order parameter is identical to the behavior found by time-dependent simulations |ll|, Il2| 
departing from a unique initial seed. In both methods it is necessary to place the system in 
its critical point. 



II. TRANSITION PROBABILITY AND SYMMETRY 

In this section and in the following we develop the formalism leading to the expression 
For this purpose we consider the class of Markovian processes defined on a lattice whose 
probability distribution P(a, t) satisfies the equation 

P(a,t)=Y / T( ( r,a',t)P (a') (4) 

cr' 

where T(a, a', t) is the probability of reaching the configuration a from configuration a' in 
an interval of time t and Pq{(t') is the initial probability distribution. We use the notation 
a = (o"i, o"2, <Jn) where N is the number of sites in the lattice and <7; is the random variable 
attached to the i-th site and that takes two values. 

If the system evolves in time according to a master equation (continuous time Markovian 
process) 

j t P(a, t) = £ a')P(a', t) - W(a', a)P(a, t)} (5) 

cr' 

then the transition probability T(a, a', t) are the elements of the matrix T given by 

f = exp{W} (6) 
where W is the matrix whose elements are 

W(a',a) = W(a',a) a' ^ a (7) 

and 

W(a,a) = - £ W(a',a) (8) 

Let R be a symmetry operation that leaves the transition probability invariant, or, equiv- 
alently, the matrix W invariant, that is 

W(Ra, Ra') = W(a, a') (9) 

and by consequence 

T(Ra, Ra', t) = T(cr, a', t) (10) 



For simplicity, we will consider here only models in which the symmetry operation R changes 
the sign of the order parameter, that is, 

M(Ra) = -M(a) (11) 

and that M(a) is linear in a, that is, 

M(<7)=£><7; (12) 

i 

III. TIME-DEPENDENT BEHAVIOR 

We will focus on the time-dependent behavior of the average 

(M(a)) t = J2M(a)P(a,t) (13) 
of the order parameter M(a). Its time evolution is given by 

(M(a))t = E E M(a)T(a, a', t)P (a') (14) 

<y a' 

where the initial state Po(c) is uncorrelated with a nonzero magnetization. That is, the 
initial magnetization 

(M(<r))o = E M (°) P M = Nm o (15) 

is nonzero, where iV is the number of sites of the lattice and m is a small quantity As 
stated by the short-time scaling theory, the order parameter folows, at the critical point, a 
power law behavior 

(M(<j)) t ~ m t e (16) 

for small values of m . According to this theory yet the initial state is prepared in such a 
way that all sites are uncorrelated with a nonzero (and small) initial order parameter mo- In 
order to set up such an initial state, one attributes to each site a magnetization mi = mofJ>i- 
Or equivalentely, the spin of the i-th site will be <7j = Hi with probability (l + m )/2 and will 
be <7j = —fj, with probability (1 — m )/2. The initial probability Po(c) can then be written 
as 

Po(cr) = + mow} (17) 

i 

where 

= ^ (18) 



Notice that using equation ()12j) and ([17]) we can trivially find that 

(M(a)) = E M(a)P (a) = Nam (19) 

where a is the constant 

a = 

N 

For small values of mo, the expansion of the initial probability Pq{<j) in powers of mo 



(20) 



gives, up to linear terms in mo the following expression 

P (a) = $ {1 + m M(a)} (21) 
Substituting this expression in equation (JT3j) we get 

(M(a)) t = E E M(a)T(a, a', t)%+ 

a a' 

+ E E M(a)T(a, a\ t)$ m M(a') (22) 

a a 1 

Now, the first term on the right hand side vanishes identically due to the following 
reasoning. Since the symmetry operation R leaves the transition probability invariant but 
changes the sign of the order parameter, we have 

E E M(Ra)T(Ra, Ra 1 , t)$ = - E E M(a)T(a, a', t)% (23) 

c a' a a' 

By a change of variable, Ra — > a, the left-hand side of this equation equals the first term of 
the right-hand side of equation (|22j) so that it vanishes. Therefore 

Q(t) = Um iMMl = ^2Y: M(a)T(a, a', t)M(a')® (24) 
and, from equation (fTBj) it follows that 

Q(t) ~ t e (25) 

IV. APPLICATIONS 

A. Models with up-down symmetry 

We begin with a simple example, namely the ferromagnetic Ising model coupled to a 
stochastic dynamics such as the Metropolis algorithm. The order parameter is defined by 

M(<7) = 5>< (26) 
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where the summation is over all sites of the lattice. For the present case fii = +1 for all 
sites of the lattice. Here the symmetry operation R, with the properties given by equations 
(|1(J|) and (J lip , is the one in which the up-down symmetry is observed, that is, the operation 
that changes the sign of each spin variable <7j — ► —a*. 

The short-time behavior for the Ising model has been already studied through the present 
formalism ^ . Besides, using this formalism, it has been possible to determine the short-time 
behavior for several irreversible models (i. e., lacking detailed balance) ^, 14] with up- 



down symmetry dynamics. These include for instance the majority vote model and similar 
nonequilibrium models jsl 



B. Models with antiferromagnet ordering 

In this case the system is divided into two sublattices A and B. The order parameter is 
defined by 

M(a) =5>i-$>< (27) 

Therefore, for this case one has //j = +1 if i £ A and \ii — — 1 if i £ B. The symmetry 
operation R is a translation such that a given site of one sublattice goes into a site of the 
other sublattice. 



C. Baxter-Wu model 



We consider in this subsection the Baxter-Wu model with triplet interactions defined on 



a triangular lattice |9|,ll0|. The lattice is composed of three sublattices which we denote by 
A, B, and C. The Baxter-Wu model does not have a global symmetry but semi-global 
symmetries. The Hamiltonian of the model and a fortiori the transition probability is 
invariant if we change the signs of two sublattices leaving the third invariant. It is convenient 
therefore to define three symmetry operations, denoted by Ra, Rb and Rc- The symmetry 
operation Ra changes the signs of the spins belonging to the sublattices B and C and leaves 
the signs of the spins of sublattice A invariant. Similar definitions can be stated for Rb and 
Rc- Each of these symmetry operations leaves the Baxter-Wu Hamiltonian invariant and a 
fortiori the transition probability invariant. 
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We take as the order parameter the magnetization of one of the sublattices, say, sublattice 
A, given by 

M A (<T) = y £a i (28) 

Comparing it with equation (fT2"|) we have that jii — 1 if i £ A and Hi = if i G B or i 6 C. 
The symmetry operation Rb (or Re) changes the sign of M A {(?) and leaves the transition 
probability invariant. According to the formalism developed in the previous section we 
conclude that 

Q A (t) = lim (MA(a))t = E E ^(a)T(a, a', i)M^(a')$o (29) 

m a J 

will behave as 

~ t 6 (30) 

Equivalently, we may demonstrate that the analogous quantities Qs(t) and Qcif) related 
to the magnetizations Mb(c) and Mc(<r) of sublattices B and Re, respectively, will behave 
as t e . 

We may also use as the order parameter the total magnetization 

M(<t) = 5>< (31) 

i 

which we write as the sum of the magnetizations of the three sublattices 

M(a) = M A {a) + M B {a) + M c (a) (32) 

which leads to 

Q(t) = E E ^)2>, ^, t)M(a')% (33) 

Substituting ((H2|) into (jHEJ) we see that is a sum of nine terms. The terms that involve 
magnetizations of distinct sublattices will vanish. For instance, the term that involve M A 
and Mb will change sign by the use o the symetry operation R A . The nonvanishing terms 
are the three terms that involve the same magnetization. One concludes that the quantity 

Q(t)=Q A (t) + Q B (t) + Qcit) (34) 

and therefore will behave as t e . This procedure was used to determine the exponent 
6. The numerical results give very precise values for the exponent when compared with the 
results coming from simulations with nonzero initial magnetization. 



V. CONTACT MODEL 



Now we discuss the short-time behavior fo models with an absorbing state. These models 

do not possess symmetry operations like the ones defined in the preceeding sections. Due 

to this fundamental difference we need to proceed by introducing another approach. The 

I I 

simplest example of this type of model is the contact process [12j. Such model is defined in 
a lattice and each microscopic state is identified with a = (a±, cr 2 , crjy) where cr, = or 1 
according wether the site i is empty or occupied by a particle. It evolves in time according to 
local Markovian rules where particles are catalitically created and spontaneously anihilated. 

The initial probability is such that all sites are uncorrelated and such that the average 
(o-j) = po, that is, 

Po{<?) = TIH 1 - <*i) + ben} (35) 

i 

where 

a = 1 - po b = po (36) 

is the total number of sites in the lattice. 

Following the short-time scaling theory the order parameter (11(a)) given by 

(n(a)) t = J2J2n(a)T(a,a',t)P (a') (37) 

a 1 a 

where 

n(<r)=£> (38) 

i 

is the number of particle, behaves, at the critical point, as 

(n(a)) ~ Po t e (39) 

Consequently the quantity 

Q( t) = I Ito (40) 
has a similar behavior in the early time regime 

Q(t) ~ t e (41) 

For small values po of we have 

P (a) = ^^o^ + a^b^^Acr) (42) 



where 

M<r) = UQ--<rt) (43) 

i 

is the probability distribution such that the configuration (0, 0, 0, 0) (all sites empty) has 
probability 1 and the other configurations have zero probability, 

*» = *i IK 1 -*) ( 44 ) 

is the probability distribution such that the configuration (0, 0, 1, ...0) (a particle placed 
at the j'-th site and all other sites empty) has probability 1 and all other configurations have 
zero probability. 

The average (n(a))t can then be written as a sum of two parts 

("(*)>* = a N E E ^ *)*o(<0 + 

+ ^ EEE "(^, a', t)$, (a') (45) 

3 cr' CT 

Since the contact process has an abosrbing state devoided of particles, the first term vanishes 
identically because $o(c) is the absorbing state. Therefore, using the translational invariance 
we obtain 

Q(*) = EE n Wn^.W ( 46 ) 

cr' cr 

To calculate numerically Q(t), we start from a configuration with just one occupied site and 
determine the number of occupied sites at time t. 

According to the scaling relations established for the time-dependent behavior of the 
contact model in which the simulation is started with just one occupied site, the average 



number of particles n p (t) behaves as 
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n p (t) ~ f (47) 

As Q(t) is identified with n p (t) so the exponent 9 is identified with the exponent r\ [16]. 
Let us now calculate the time correlation of a given site, say the site j. It is given by 

A(i)=EE#. fl '.«)W ( 48 ) 

cr' cr 

and behaves, according to the short-time scaling theory, as 

A(t) ~ t x (49) 



Now 

(jjPoia) = bo 3 [] " <?i) + 6^0 (50) 

where we have used the obvious relations 0,(1 — Cj) = and Oj<Jj = o r Therefore, in the 
limit po - *• 0, we get 

lim ^P (a) = $ J (a) (51) 



Consequently, 



B(f) = hm ^ = E E ^(a, a', t)$, (a') (52) 



so that -B(t) behaves as 

B(t) ~ £ A (53) 

Given that the initial particle seed is placed at a given site, the quantity B(t) is the proba- 
bility that this site be occupied at time t. 

The exponent A is related to dynamic exponent z by A = d/z — 9 [8]. Since the exponent 
9 was identified as the exponent rj, it follows that A = d/z — p. Now, from the hyperscaling 
relation for the contact process we have d/z — r/ = 25 [12j where 5 is the exponent associated 
to the survival probability. Therefore it follows that A = 25. 

^From the formalism developed here we conclude that the study of the short-time behavior 
of the contact process, discussed in reference [l(| , is as a matter of fact equivalent to the study 
of this model by means of the time- dependent technique. Moreover, the critical exponents 
associated to the short time dynamics for the contact model, as well as the relation among 
them, are equivalent to those found for the time-dependent simulations. 

VI. CONCLUSION 

We have shown that the short-time critical exponent 9 of several models invariant under a 
given group of symmetry can be calculated from the time correlation of the order parameter 

Q{t) = (M(t)M(O)) (54) 

where (f(t)g(0)) is a notation defined by 

(f(t) g m = E E n°)n^ t)gw)*» (55) 
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We have also obtained similar formula for the contact process and shown that the short- 
time critical exponent 9 is equal to time-dependent critical exponent rj. Finally, the results 
obtained here for the continuous time Markovian processes can be straitforwardly extended 
to the probabilistic cellular automata (discrete time Markovian process). 
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